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ABSTRACT 

Intrinsic alignment (IA) of source galaxies is one of the major astrophysical systematics for 
ongoing and future weak lensing surveys. Several IA models have been proposed in the litera¬ 
ture and constrained through observations. Their impact on cosmological constraints has been 
examined using the Fisher Information Matrix in conjunction with approximate covariance 
computation schemes. 

This paper presents the first forecasts of the impact of IA on cosmic shear measurements for 
future surveys (DES, Euclid, LSST, WFIRST) using simulated likelihood analyses and realis¬ 
tic covariances that include higher-order moments of the density field in the computation. We 
consider a range of possible IA scenarios and test mitigation schemes, which parameterize IA 
by the fraction of red galaxies, normalization, luminosity and redshift dependence of the IA 
signal (for a subset we consider joint IA and photo-z uncertainties). 

Compared to previous studies we find smaller biases in time-dependent dark energy models if 
IA is ignored in the analysis; the amplitude and significance of these biases vary as a function 
of survey properties (depth, statistical uncertainties), luminosity function, and IA scenario: 
Due to its small statistical errors and relatively shallow observing strategy Euclid is signifi¬ 
cantly impacted by IA. LSST and WFIRST benefit from their increased survey depth, while 
the larger statistical errors for DES decrease IA’s relative impact on cosmological parameters. 
The proposed IA mitigation scheme removes parameter biases due to IA for DES, LSST, and 
WFIRST even if the shape of the IA power spectrum is only poorly known; successful IA 
mitigation for Euclid requires more prior information. We explore several alternative IA mit¬ 
igation strategies for Euclid; in the absence of alignment of blue galaxies we recommend the 
exclusion of red (IA contaminated) galaxies in cosmic shear analyses. We find that even a re¬ 
duction of 20% in the number density of galaxies only leads to a 4-10% loss in cosmological 
constraining power. 
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1 INTRODUCTION 

Weak gravitational lensing by large-scale structure, so-called cos¬ 
mic shear, is a promising technique to constrain cosmological pa¬ 
rameters. Cosmic shear is among the key science cases of vari¬ 
ous ongoing surveys, such as Kilo-Degree Survey (KiDS 1 ), Hyper 
Suprime Cam (HSC 2 ), and Dark Energy Survey (DES 3 ). Indepen¬ 
dent of any assumptions about the relationship between dark and lu¬ 
minous matter, it provides valuable information on both the geome- 

* E-mail: aelisabeth.krause@gmail.com 

1 http://www.astro-wise.org/projects/KIDS/ 

2 http://www.naoj.org/Projects/HSC/HSCProject.html 

3 www.darkenergysurvey.org/ 


try and structure growth of the Universe (Hoekstra, Yee & Gladders 
2002; van Waerbeke, Mellier & Hoekstra 2005; Jarvis et al. 2006; 
Schrabback et al. 2010; Lin et al. 2012; Heymans et al. 2012; Huff 
et al. 2014). 

If systematics can be sufficiently controlled for future mis¬ 
sions, such as the Large Synoptic Survey Telescope (LSST 4 ), Eu¬ 
clid 5 and the Wide-Field Infrared Survey Telescope (WFIRST 6 ), 
cosmic shear has the potential to be the most constraining cos¬ 
mological probe. If systematics control does not improve, lensing 

4 http://www.lsst.org/lsst 

5 sci.esa.int/euclid/ 

6 http://WFIRST.gsfc.nasa.gov/ 
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constraints will be substantially weaker than constraints from other 
probes, such as SNla, Baryonic Acoustic Oscillations, and redshift- 
space distortions (Weinberg et al. 2013). 

Important systematic effects that complicate the extraction of 
cosmological information from cosmic shear are shear calibration 
(Hirata & Seljak 2003; Huterer, Keeton & Ma 2005) and photo-z 
errors (Ma, Hu & Huterer 2006; Bernstein & Huterer 2010; Hearin 
et al. 2010), baryonic effects (Jing et al. 2006; Zentner, Rudd & 
Hu 2008; Semboloni et al. 2011; Zentner et al. 2013; Semboloni, 
Hoekstra & Schaye 2013; Eifler et al. 2014a) and intrinsic align¬ 
ment (IA) of source galaxies. 

Cosmic shear is typically measured through two-point corre¬ 
lations of observed galaxy ellipticities. In the weak lensing regime, 
the observed ellipticity of a galaxy is the sum of intrinsic ellip- 
ticity e 1 and gravitational shear y: e obs » e l + y. If the intrinsic 
shapes of galaxies are not random, but spatially correlated, these 
intrinsic alignment correlations can contaminate the gravitational 
shear signal and lead to biased measurements if not properly re¬ 
moved or modeled. Since early work establishing its potential ef¬ 
fects (Croft & Metzler 2000; Heavens, Refregier & Heymans 2000; 
Catelan, Kamionkowski & Blandford 2001; Crittenden et al. 2001), 
IA has been examined through observations (e.g., Hirata et al. 
2007; Joachimi et al. 2011; Blazek et al. 2012; Singh, Mandelbaum 
& More 2014), analytic modeling, and simulations (e.g., Schneider, 
Frenk & Cole 2012; Tenneti et al. 2014b,a) - see Troxel & Ishak 
(2014), Joachimi et al. (2015), and references therein, for recent re¬ 
views. A fully predictive model of IA would include the complex 
processes involved in the formation and evolution of galaxies and 
their dark matter halos, as well as how these processes couple to 
large-scale environment. In the absence of such knowledge, ana¬ 
lytic modeling of IA on large scales relates observed galaxy shapes 
to the gravitational tidal field, and typically considers either tidal 
(linear) alignments, or tidal torquing models. 

The shapes of elliptical, pressure supported galaxies are of¬ 
ten assumed to align with the surrounding dark matter halos, which 
are themselves aligned with the stretching axis of the large-scale 
tidal field (Catelan, Kamionkowski & Blandford 2001; Hirata & 
Seljak 2004). This tidal alignment model leads to shape alignments 
that scale linearly with fluctuations in the tidal field, and it is thus 
sometimes referred to as “linear alignment,” although nonlinear 
contributions may still be included (Bridle & King 2007; Blazek, 
McQuinn & Seljak 2011; Blazek, Vlah & Seljak 2015). For spi¬ 
ral galaxies, where angular momentum is likely the primary factor 
in determining galaxy orientation, IA modeling is typically based 
on tidal torquing theory, leading to a quadratic dependence on tidal 
field fluctuations (Pen, Lee & Seljak 2000; Catelan, Kamionkowski 
& Blandford 2001; Hui & Zhang 2002; Lee & Pen 2008), although 
on sufficiently large scales, a contribution linear in the tidal field 
may dominate. Due to this qualitative difference in assumed align¬ 
ment mechanism, source galaxies are often split by color into “red” 
and “blue” samples, as a proxy for elliptical and spiral types. In¬ 
deed, blue samples consistently exhibit weaker IA on large scales, 
supporting the theory that tidal alignment effects are less promi¬ 
nent in spirals (Faltenbacher et al. 2009; Hirata et al. 2007; Man¬ 
delbaum et al. 2011). On smaller scales, IA modeling must in¬ 
clude a one-halo component to describe how central and satellite 
galaxies align with each other and with respect to the distribution 
of dark matter (Schneider & Bridle 2010). Numerical simulations, 
especially those including hydrodynamical physics, have recently 
become powerful tools for constructing these models (Schneider, 
Frenk & Cole 2012; Joachimi et al. 2013b,a; Tenneti et al. 2014b,a). 
However, due to the cumulative uncertainty in modeling the cosmic 


shear signal on angular scales corresponding to the one-halo term 
regime, e.g. from uncertainties in modeling the non-linear matter 
power spectrum (e.g., Heitmann et al. 2009; Eifler 2011) and the 
impact of baryons (e.g., van Daalen et al. 2011; Zentner et al. 2013), 
a more conservative analysis strategy might be to exclude corre¬ 
sponding scales from the analysis (e.g. Faureijs et al. 2011). 

In this work, we focus on IA for red source galaxies on scales 
outside the one-halo regime, testing the potential impact on cos¬ 
mic shear analyses using several different tidal alignment scenar¬ 
ios. Given recent observational results, IA for blue source galaxies 
is likely subdominant. Although we examine the potential effects 
of tidal alignment for blue galaxies, we leave a detailed treatment 
of IA contamination from blue galaxies for future work. 


2 INTRINSIC ALIGNMENT MODELING 

We consider a class of shape alignment models for which the cor¬ 
related component of galaxy shapes is proportional to the gravi¬ 
tational tidal field. These models are often referred to as “linear 
alignment” or “tidal alignment,” and are most frequently used to 
describe IA of elliptical (red) galaxies. In this work, we do not use 
a separate model for the (much weaker) alignment of blue galaxies. 
For further discussion, see Catelan, Kamionkowski & Blandford 
(2001); Hirata & Seljak (2004); Joachimi et al. (2011); Blazek, Mc¬ 
Quinn & Seljak (2011); Blazek, Vlah & Seljak (2015). 

The underlying principle of tidal alignment models is that the 
gravitational collapse of an initially spherical overdensity in a tidal 
gravitational field leads to triaxial haloes, such that the halo will be 
prolate if the overdensity is stretched by the large scale tidal field 
and oblate if it is compressed (Catelan, Kamionkowski & Bland¬ 
ford 2001). Under the assumption that the shape and orientation 
of an elliptical galaxy are determined by the shape of the halo in 
which it resides, this mechanism leads to a net correlation of intrin¬ 
sic galaxy ellipticities with the gravitational tidal field. The lowest- 
order contribution is a linear relationship to the tidal field, which 
can be expressed in Fourier space in terms of the density contrast 
5: 

(k : lc*,2k x ky) 

e^(k,z) = -A(z)— - k2 S [<5(k,z IA )j , (1) 

where the ellipticity components, denoted (+, x) are measured with 
respect to the x-axis. The amplitude A(z ) depends on redshift and 
galaxy properties, and S is a filter that smooths the tidal field on 
the relevant scale (see Blazek, Vlah & Seljak 2015). The redshift 
at which IA is set, zia, the amplitude A(z), and the treatment of the 
density field can differ in various implementations of the tidal align¬ 
ment model, as detailed below. The intrinsic alignments contribute 
to the observed ellipticity correlations through through both their 
auto-correlation and their cross-correlation with the density field, 
known respectively as the “II” and “GI” terms (see Hirata & Sel¬ 
jak 2004). At leading order, gravitational lensing produces £-mode 
(curl-free) shape correlations, and we are thus concerned with the 
£-mode IA contribution. 7 For simplicity, we introduce the follow¬ 
ing notation for the E-mode component of the intrinsic ellipticity 
field: 

4 A = -A(z)£(Lz). (2) 

7 5-modes can be produced by IA contamination as well as higher-order 
lensing effects. 
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The relevant 3D IA power spectra can then be written: 


2.2 IA density power spectra models 


Pn(k,z)=A\z)P SB (.k,z), (3) 

P G i(k,z) = -A(z)Pss(k,z). (4) 


2.1 IA amplitude 


The redshift at which the alignment is set, zia, and the subse¬ 
quent IA evolution are determined by the astrophysical processes 
involved in galaxy formation and evolution. It is sometimes as¬ 
sumed that zia is during matter domination when a halo first forms 
(“primordial alignment”). However, late-time accretion and merg¬ 
ers could have a significant impact on IA, in which case the relevant 
Zia could be closer to the observed redshift, z 0 b s . In the case of in¬ 
stantaneous alignment (zia = Zobs), 

A tast (z) = C“ st pm,o(l+z). (5) 


with a proportionality constant Cj describing the strength of galaxy 
shape alignments, which may depend of galaxy properties and 
needs to be determined empirically (see Blazek, Vlah & Seljak 
2015, for further discussion on A(z)). 

In the “primordial alignment” assumption, IA is set at some 
early time and does not evolve. In this case, since the relevant den¬ 
sity field is determined at a redshift where nonlinear growth can be 
neglected, we can instead evaluate it at the observed redshift and 
absorb the linear growth factor, D(z) into the overall amplitude: 

A early (z) = OAo(l + z,a)^ = C ‘ , (6) 

D(z) D(z ) 

where in the second step we have assumed that alignment happens 
during matter domination such that M = (1 + zia)D(zia) * s constant. 
The resulting expression has no dependence on zia- 

In practice, the redshift evolution of intrinsic alignments may 
be between these scenarios, and the effective amplitude will depend 
on the full luminosity distribution of the (red) galaxy sample of 
interest. In this analysis, we use the redshift and luminosity scaling 


A(L,z) = 


ClPm.O 

D(z) 


Ao 




(V) 


as baseline model. In this expression Cjp cr ~ 0.0134 is a normal¬ 
ization derived from SuperCOSMOS observations (Hirata & Seljak 
2004; Bridle & King 2007), and for the fit parameters A, t/ 0 thcr,j 8 
we adopt the constraints from the MegaZ-LRG + SDSS LRG sam¬ 
ple in Joachimi et al. (2011), i.e. A 0 = 5.92 y °^, 77 = —0.47^^, 
and [3 = 1. IO^q Iq, with zo = 0.3 the observationally-motivated 
pivot redshift and L 0 the pivot luminosity corresponding to an ab¬ 
solute r-band magnitude of -22. For our implementation, we aver¬ 
age A(L, z) over the red galaxy luminosity distribution to compute 
the average IA amplitude for the red source galaxy population at 
a given redshift (see Sect. 4.1), and include an additional parame¬ 
ter 77high-z t° account for the uncertainty in the extrapolation of this 
redshift scaling, 


©(zi - z) + 0(z - Zi) 

( 8 ) 

with m Um the survey’s limiting magnitude, and 0 the step function. 
The term in square brackets is a truncated power-law in (1 + z), 
which is unity for z <= zi and with slope phigh-z for z > zi- As 
the fiducial redshift scaling is based on the MegaZ-LRG + SDSS 
LRG sample with z < 0.7, we choose zi = 0.7 and marginalize over 
^high-z (see Sect. 5.1 and Tab. 4 for details). 


A(m lim , z) = (A(L, z))^ ^ x 


/ 1 + z \"“ 8h - z 

' 1+zJ 


Four effects can produce nonlinearities in the intrinsic shape corre¬ 
lations: ( 1 ) nonlinear dependence of intrinsic galaxy shape on the 
tidal field (e.g., quadratic tidal torquing); ( 2 ) nonlinear evolution 
of the dark matter density field, leading to nonlinear evolution in 
the tidal field; (3) a nonlinear bias relationship between the galaxy 
and dark matter density fields; (4) the IA field actually observed is 
weighted by the local density of galaxies used to trace the shapes. 
In this work, we assume a pure tidal alignment model, i.e. that in¬ 
trinsic galaxy shapes depend linearly on the tidal field, even on 
small scales. However, nonlinearities from the other three effects 
must still be considered. The following versions of the tidal align¬ 
ment model treat these effects differently. We also implement tidal 
field smoothing only in the full tidal alignment model, as described 
below, allowing us to test the impact of different smoothing proce¬ 
dures, since the physically correct smoothing is unknown. 

Linear Alignment (LA) The most basic version of tidal alignment 
includes only the leading order contribution to IA correlations (i.e. 
from the linear density field dn n ): 

Psdk,z) = P«:(Lz) = P\m(k, z) (9) 

Nonlinear Alignment (NLA) To account for the nonlinear evolu¬ 
tion of the density field, the linear power spectrum in Eq. 9 can be 
replaced with the nonlinear version: 

Pss(k, z) = Pss(k, z) = Pudk, z) (10) 

This approach, often called the nonlinear alignment (NLA) model, 
was first used by Hirata & Seljak (2004) and more fully developed 
by Bridle & King (2007). Due to its relative simplicity and abil¬ 
ity to more accurately model observed intrinsic shape correlations, 
the NLA model has been used in weak lensing forecasts and analy¬ 
ses (e.g. Heymans et al. 2012). However, the NLA model neglects 
nonlinearities from source density weighing and nonlinear galaxy 
biasing, the former of which can be similar in size or larger than 
the correction due to nonlinear structure growth. 

Freeze-in Model (FR) In the primordial alignment scenario, in¬ 
trinsic shape correlations are set at early times, when nonlinear 
evolution of the density field can be largely ignored, and are thus 
proportional to the linear tidal field. For the GI term, however, the 
nonlinear growth of the density field should still be included. To 
capture these effects, Kirk et al. (2012) suggested a model in which 
the density-ellipticity cross-correlation was given by the geometric 
mean of the linear and nonlinear power spectra: 

Pss(k,z ) = dPnn(k,z)PNL(k,z), (11) 

Pes(k,z) = PyJk,z). (12) 

Although approximate, this approach allows the cross-correlation 
of the density field at different redshifts, without needing to account 
for the effects of advection that appear in an Eulerian perturbation 
theory approach (see (Blazek, Vlah & Seljak 2015) for further dis¬ 
cussion). Note that we normalize all IA models to match the ob¬ 
served large-scale Rgi amplitude at the pivot redshift zo = 0.3 of 
Joachimi et al. (2011), hence our A(z) differs from the normaliza¬ 
tion chosen in Kirk et al. (2012). After accounting for the correct 
redshift dependence of A(z), our implementation of the freeze-in 
P n model exhibits same redshift scaling as Eq. ( 8 ) in Kirk et al. 
(2012), and our P GI expression differs by one power of the linear 
growth function, which is missing in their Eq. (9). 
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Full tidal alignment model (TA) Although the NLA approach im¬ 
proves the model fit to data, it is not fully consistent and omits 
important astrophysical effects. The model of Blazek, Vlah & Sel- 
jak (2015) includes all terms that contribute at next-to-leading or¬ 
der while simultaneously smoothing the tidal field (e.g. at the La- 
grangian radius of the host halo). The effects of weighting by the 
local galaxy density can be larger than the correction from the non¬ 
linear evolution of dark matter density, especially in the case of 
highly biased tracers, since the leading effect scales with the linear 
galaxy bias. This more complete tidal alignment model is given by 
(Blazek, Vlah & Seljak 2015): 

58 , 

Pss(k,z) - Pnl(.Kz) +-^^bicr s P\j D +biPo\os, (13) 

116 , ,, 

Psa(k,z) - T’nlC^jZ) + + 2fcl/ > 0|06 + ^l^ , 06|06>(14) 

where b\ is the linear bias of the source sample, and Poioe and /’oeioe 
are 0(P 2 in ) terms that arise from weighting the intrinsic shape field 
by the source density. In these expressions, a Gaussian smoothing 
filter with scale length k sm = 1.0/7~'Mpc is applied to the tidal field 
and consistently treated in the subsequent calculations. The density 
variance cr| is determined by integrating the linear power spectrum 
smoothed with this filter. For this IA scenario, the redshift and lu¬ 
minosity scalings of the amplitude are determined by the evolution 
of b i and cr^, rather than applying Eq. (7). The fiducial value of Ao 
at pivot redshift zo is multiplied by (1 + 58/105ii(zo)c|(zo)) _1 to 
ensure that P GI has the correct amplitude on large scales at that red¬ 
shift (where measurements provide the tightest constraints). Note 
that P\\ will not have the same large-scale amplitude as in the other 
models, since the density-weighting correction is different. 


3 MODELING WEAK LENSING OBSERVABLES AND 
COVARIANCES 

All simulated likelihood analyses in this paper are computed using 
the weak lensing modules of CosmoLike (see Eifler et al. 2014b, for 
an early version; official release paper Krause et al. 2015 in prep). 


Wo + (1 - a)w a ) are incorporated following the recipe of icosmo 
(Refregier et al. 2011), which in the non-linear regime interpo¬ 
lates Halofit between flat and open cosmological models (also see 
Schrabback et al. 2010, for more details). 

Having obtained the density power spectra we calculate the 
shear power spectra as 


C\l) : 


qw 4 o 2 
4 C 4 


f 


g(X)g J <X) n 
<k — T-— p s 

a 2 (x) 


f K (xY 


x 


(15) 


with / being the 2D wave vector perpendicular to the line of sight, 
X denoting the comoving distance, Xh is the comoving distance to 
the horizon, a(x) is the scale factor, and f' K (x ) the comoving angu¬ 
lar diameter distance (throughout set to x since we assume a flat 
Universe). The lens efficiency g‘ is defined as an integral over the 
redshift distribution of source galaxies n(x(z)) (see Sect. 4 for de¬ 
tails) in the i th tomographic interval 

r x '' A, f v' — y) 

g'(T)= Ax'n'ix' Y * ■ G6) 

We split the source galaxies into five tomographic bins, chosen such 
that each bin contains the same number of galaxies. 


3.2 Projected intrinsic alignment power spectra 

For a given a intrinsic alignment model, we compute the corre¬ 
sponding projected power spectra as 

^ ij n , f Xh 4 M'CfVUO j-2 i , n „ , 

C n (/, rn lim ) = dx -:-/ rcd ( ,n ii m , z) P\i (k , z, rniim) , 

Jo X 

(17) 

n ... C Xh 4 .yCr)g%) r _ n ...... s 

C Gl (l, m\im) ~ I hV 2 Ae d(^lim? w \A> ^llim) ■> 

Jo X 

(18) 

where / re d, the fraction of red galaxies at redshift z(x ) in a particular 
survey, is computed from the luminosity functions for all and red 
galaxies, as described in Sect. 4.1. 


3.1 Shear tomography power spectra 

We compute the linear power spectrum using the Eisenstein & 
Hu (1999) transfer function and model the non-linear evolution of 
the density field using the updated Halofit as described in Taka- 
hashi et al. (2012). Time-dependent dark energy models (w = 

I 


3.3 Shear covariances 

The non-Gaussian covariance of tomographic shear power spectra 
C‘J(l) can be written as (Cooray & Hu 2001; Hu & Jain 2004; Sato 
et al. 2009) 


Cov(c' 2 (/,),Cf(/ 2 )) = 


2nd, 


1112 


njiAh 


_2 V 


c?(/o + S ik 


In 1 


r d 1 1 

r d 2 v 

J|i|£(i Mh) 

J |l']et 2 Mh) 


gal 

l 

a" 


2 nr. 


J2 \ 


C''(/,) + 6 a 


2 n 


gal / 


Ci k (h) + 6 jk —f 
V 2W gaD 


r^(i,-u',-T) + r^s V (i,-i,i',-T) 

I- 


(19) 


with fl s the survey area, n'. the number of source galaxies in to¬ 
mography bin i, cr € the ellipticity dispersion, A(/;) = L d 2 1 ~ 
2nIjAlj the integration area associated with a power spectrum bin 
centered at /, and width A/,, and and T^hsv the convergence 
trispectrum in the absence of finite volume effects and the super¬ 
sample variance contribution to the trispectrum (Sato et al. 2009; 


Takada & Hu 2013). We use a halo model implementation of these 
terms as described in (Eifler et al. 2014a). 
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Figure 1 . Normalized redshift distributions for the four different surveys. 
Top panel: DES and Euclid. Bottom panel: LSST and WFIRST. 


Table 1 . Survey parameters (see text for details). 


Survey 

£2 S [deg 2 ] 

cr e 

w gal 

£max 

2mean 

£med 

m r, lim 

DES 

5,000 

0.26 

10 

2.0 

0.69 

0.62 

24.0 

LSST 

18,000 

0.26 

26 

3.5 

1.07 

0.93 

27.0 

Euclid 

15,000 

0.26 

20 

2.5 

0.8 

0.74 

24.5 

WFIRST 

2,200 

0.26 

45 

4.0 

1.27 

1.11 

28.0 


4 SURVEY PARAMETERS 

For our IA study we consider four different surveys: DES, Euclid, 
LSST and WFIRST. The corresponding redshift distributions are 
modeled as 


n(z) = N z a exp 



( 20 ) 


The redshift distributions are shown in Fig. 1 with mean and me¬ 
dian redshift summarized in Table 1. For LSST we adopt the red¬ 
shift distribution suggested in Chang et al. (2013), and the DES 
redshift distribution is modeled by a modified CFHTLS redshift 
distribution (see Benjamin et al. 2007, adjusted for the slightly 
lower mean redshift of DES). The Euclid redshift distribution fol¬ 
lows a similar parameterization as DES, but is assumed to be 
slightly deeper. WFIRST redshifts are modeled after the LSST z- 
distribution, again assuming slightly deeper imaging. It is impor¬ 
tant to note that throughout this paper WFIRST refers to the wide- 
held imaging component of the WFIRST High-Latitude Survey 
(HLS) as described in Spergel et al. (2015). The HLS component 
of WFIRST encompasses only 2 years and includes a substantial 
amount of spectroscopic grism observations. 


Table 2. Fiducial luminosity function parameters. 



^[(/7/Mpc) J ]f 

M,!(0.1) t 

a r 

P T 

e + 

pIM 

e h * 

all 

9.4 X 10 -3 

-20.70 

-1.23 

1.8 

0.7 

-0.30 s 

1.23 

red 

1.1 x 10~ 2 

-20.34 

-0.57 

-1.2 

1.8 

—1.15 s 

1.20 


f () l r-band LF fit parameters from Tab. 5 in Loveday et al. (2012). 
f B-band LF evolution parameters from Tab. 6 in Faber et al. (2007). 
s Value rescaled to match definition in Eq. (22). 


4.1 Source galaxy distribution 


We use an evolving Schechter luminosity function (LF) </>(L, z) to 
model the LF of source galaxies 

t (L ’Z ) = fv(zy ex p(-z|))’ (21) 

with evolving Schechter parameters 

<l'(z) = 4>' 0 10 04ft , (22) 

and L*(z) the luminosity corresponding to 

M*(z) = M*(0.1) - Q(z- 0.1) . (23) 


For the baseline model this analysis uses the r-band fit parameters 
from the GAMA survey (Loveday et al. 2012), which are listed in 
Tab. 2. 

As described in Sect. 2.1, we follow Joachimi et al. (2011) in 
modeling the luminosity and redshift dependence of the IA ampli¬ 
tude. We evaluate the effective amplitude (A) and the fraction of 
contaminated (red) galaxies / rcd from the LF of all and red galaxies 


A(miim, z ) 


,/redOniim. Z) 


r dLiff fa ti {L,z) 
A ( Lo , Z ) - 

fL(mn m ,z) dL< / > redL, ~) 


(24) 

(25) 


with z) the L-corrected limiting luminosity at redshift z in a 

survey with limiting apparent magnitude m u m (c.f. Tab. 1 for values 
used in this analysis to model different surveys). L(mn m ,z) is the 
luminosity corresponding to 


Aflimfc mm) = "him - |5 log 10 + 25 + L(z)j , (26) 

where ZJ L (z) is the luminosity distance, and k(z) the k+ e-correction, 
for which we use the r-band values from (Poggianti 1997) due to 
the wide redshift range of this model. 

We note that the LF parameters and k + e-correction model as¬ 
sumed in this study differ from the choices in Joachimi et al. (2011). 
While the Faber et al. (2007) LF parameter used in Joachimi et al. 
(2011) are based on deep observations LF measurements to z ~ 1, 
compared to z < 0.5 for GAMA, it is important to note that the for¬ 
mer are R-band LF parameters (compared to 01 r for the GAMA pa¬ 
rameters in our baseline model). Since Schechter parameters vary 
significantly across bands, it is not clear which of these will be 
a better model for the source galaxies of future weak lensing sur¬ 
veys with high sensitivity in the red and near-infrared bands. Hence 
the LF parameters are a major uncertainty in our forecasts for the 
unmitigated impact on IA on cosmic shear surveys. To illustrate 
the impact of LF parameters, we will consider a second LF model 
which combines the DEEP2 evolution parameters with the GAMA 
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Table 3. Fiducial cosmology, minimum and maximum of the flat prior on 
cosmological parameters. 
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0.6 
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0.6 

0.95 

1.06 

0.0 

2.5 

0.055 

0.76 


is an approximation to the correct approach of a cosmology depen¬ 
dent or estimated covariance (see Eifler, Schneider & Hartlap 2009, 
for further details). 

The posterior probability for a point in cosmological parame¬ 
ters p c is obtained via Bayes’ theorem 


Figure 2. Fraction of red galaxies computed from GAMA and DEEP2 lu¬ 
minosity function, respectively. We consider a deep and a shallow survey 
with limiting magnitude of 27.5 (LSST/WFIRST) and 24.5 (Euclid/DES), 
respectively. 


LF at the pivot redshift ( P D , Q D in Tab. 2), while this differs from 
the Faber et al. (2007) LF we will refer to as DEEP2 LF in the fol¬ 
lowing for simplicity. Figure 2 shows / re d derived from the GAMA 
and DEEP2 LFs as a function of redshift. The large discrepancy 
between the curves indicates that the LF modeling of upcoming 
and future surveys limit is a major uncertainty for current forecasts. 
Compared to the uncertainty in LF evolution, the spread between 
k + e-correction templates is a minor source of uncertainty. 


5 RESULTS 

5.1 Simulated Likelihood analyses 

Throughout this study, we quantity the impact of IA on cosmologi¬ 
cal parameters using simulated likelihood analyses as described in 
this section. We generate a IA contaminated data vector as 

D = C ij (l) + C[{(1) + C&(Z) (27) 

where the individual terms are calculated from Eqs. (15), (17), and 
(18). We split the source galaxies into five tomographic bins with 
equal number density of galaxies. From these five tomographic bins 
we derive five auto and ten cross power spectra, which we sample 
with 12 logarithmically spaced bins in / ranging from / m i n = 100 to 
/max = 5000. 

The model vector is computed in a similar way; depending 
on whether we want to examine the impact of an IA scenario or 
mitigate IA, we calculate M = CHI) or M = C ij (l) + C','(/) + Cj/,(/), 
respectively. While we contaminate the data vector with range of 
IA models (c.f. Sect. 2.2), the mitigation strategy used throughout 
this analysis only uses the NLA Halofit IA model to compute the 
IA contribution to the model data vector M. 

Using the covariance defined in Sect. 3.3 we compute the like¬ 
lihood of the (IA contaminated) data vector D given a set of cos¬ 
mological parameters p c as 

L(D|p c ) = N x exp|-i((D-M)'C _1 (D-M))|, (28) 

where we have assumed that the errors of the data vector are dis¬ 
tributed as a Multivariate Gaussian. C is the covariance matrix, 
which we approximate to be independent of any of the parame¬ 
ters. We note that assuming a constant, known covariance matrix C 


P(p c |D)ocP,.(p c )L(D|p c ), 


(29) 


where P, (p c ) denotes the prior probability. 

As discussed in Sect. 4.1, the luminosity function evolution 
parameters and the faint end slope of the luminosity function are 
key uncertainties for the amplitude of the IA contamination, and 
we treat these parameters as part of our IA model. In absence of 
reliable LF parameter constraint estimates for future surveys, we 
include abundance information to provide an implicit prior on LF 
parameters. Specifically, we include the number density of all and 
red galaxies at the mean redshift of each tomography bin, assum¬ 
ing that these will be measured with 5% uncertainty. In addition to 
this abundance prior, we reject LF parameter combinations which 
cause unphysical models with / re d > 1, which are an artifact of the 
phenomenological choice of parameterization. 

We run simulated likelihood analyses in seven-dimensional 
parameter space (see Table 3) and integrate over ten IA nuisance 
parameters (see Table 4) when the CosmoLike IA mitigation mod¬ 
ule is used. This marginalized likelihood reads 


L(D|p c )= J dp n exp^-i(D - M(p c , pJ/CT^D - M(p c , p n ))j 

(hj(z/) - rijizi, p)) 


x exp 


: I S- 

je{al!,red) i=l 


c 2 (w/fo)) 


(30) 


with p n spanning the ten nuisance parameter dimensions, and 
where hj(z.i) is fiducial galaxy abundance, nj(zi, p) the galaxy abun¬ 
dance prediction for model parameters p. 


5.2 LSST baseline models 

In this section we present detailed likelihood analyses for LSST 
assuming the NLA Halofit scenario. 

Impact of luminosity functions As detailed in Sect. 4.1, the ex¬ 
trapolation of existing LF measurements to the depth of future 
cosmic shear surveys is a major uncertainty in modeling IA. In 
Fig. 3 we illustrate the range of this uncertainty by running simu¬ 
lated likelihood analysis using the GAMA (red/dashed) and DEEP2 
(green/long-dashed) luminosity functions in the IA contaminated 
data vector. Both analyses assume the NLA Halofit scenario in the 
data vector; the corresponding contours mimic likelihood analyses 
which do not account for the IA contamination at all. As a result, 
corresponding constraints are biased at the level of several a (the 
black “x” indicates the input cosmology). The LSST statistical er¬ 
rors are shown in black/solid. We also run likelihood analyses that 
include the CosmoLike IA mitigation module, but defer the discus¬ 
sion of these results to Sect. 5.3. 
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Table 4. Left: Fiducial, minimum and maximum values (flat priors) for the intrinsic alignment parameters. Right: Fiducial values and range of the Gaussian 
priors for the nuisance parameters describing photo-z uncertainties (optimistic and pessimistic scenario) and uncertainties in the luminosity function. See text 
for details. 
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Figure 3. The impact of IA on WL constraints (68% confidence region) from LSST assuming the NLA Halofit scenario. We consider different luminosity 
functions i.e. GAMA (red/dashed) and DEEP2 (green/long-dashed) and for the GAMA LF we also consider the case for which blue galaxies have a mild NLA 
IA contribution (blue/dot-dashed). The LSST statistical errors are shown in black/solid. Orange/dot-long-dashed contours show results when using the most 
extreme of these cases, i.e. the data vector corresponding to the blue contours, as input and including the CosmoLike IA mitigation module in the analysis. 
The marginalized likelihood is obtained by integrating over a 11 -dimensional nuisance parameter space (see text for details). 


Alignment of blue galaxies While measurements of the position- 
shape correlation for blue galaxies indicate that the NLA signal 
of blue galaxies is consistent with zero (e.g., Hirata et al. 2007; 
Mandelbaum et al. 2011), the statistical uncertainty in these non¬ 
detections is substantial and does not rule out noticeable contami¬ 
nation. Hence in Fig. 3 we also illustrate the impact of NLA con¬ 
tamination for blue galaxies with an amplitude based on the 68% 
upper limit for the NLA amplitude of WiggleZ and blue, low- 
redshift SDSS galaxies in Mandelbaum et al. (2011). Assuming the 
same redshift and luminosity scaling of the blue galaxy IA ampli¬ 
tude as in our fiducial model for red galaxies, in our notation this 
corresponds to an amplitude A(Lo,z)biue = 0.25 x A(Lo,z) re d- The 


blue/dot-dashed contours in Fig. 3 correspond to a scenario where 
both red and blue galaxies are affected by IA, and the total IA am¬ 
plitude is calculated by 

z) — ./red'^redC^fimi z) + (1 Tred)Ablue(^him) z) • (31) 

This analysis again assumes the NLA Halofit scenario with the 
GAMA luminosity function, the IA model for red galaxies de¬ 
scribed in Sect. 2, and the blue galaxies NLA model motivated 
above. Equation 24 for blue galaxies hence uses A(Lo,z)biue = 
0.25 x A(L 0 ,z) red. an d we approximate tp(L,Z) b i ue with <p(L, Z) aU . 
We note that a more realistic IA model for blue galaxies involves 
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Figure 4. Top: Marginalized WL constrains (68% confidence region) from LSST when marginalizing over Gaussian photo-z uncertainties (red/dashed) and 
joint uncertainties of photo-z’s and the fiducial IA NLA Halofit model. We assume two different levels of photo-z errors resembling optimistic and pessimistic 
LSST photo-z errors. The black/solid lines again show the LSST statistical errors for comparison. Bottom: The marginalized one-dimensional posterior 
probabilities for the 12 nuisance parameters used in the optimistic and pessimistic joint photo-z and IA analysis. 


will include a contribution from the tidal torquing mechanism; we 
defer a detailed analysis to future work (Zu et al., in prep.). 

The orange/dot-long-dashed contours in Fig. 3 show results 
when marginalizing over the IA nuisance parameters as described 
in Eq. 30. In addition to the four IA and six luminosity function 
nuisance parameters (see Table 4), in this particular analysis we 
include an additional nuisance parameter describing the fractional 
amplitude of the blue IA signal (centered around the fiducial value 
of 0.25). The bias from IA contamination ( blue/dot-dashed ) is com¬ 
pletely removed in the case of the marginalized likelihood analysis; 
the broadening of the likelihood contours is noticeable but not ex¬ 
cessive (c.f. Kirk et al. 2012, 2015, for the degradation of other IA 
mitigation techniques). 


Impact of photometric redshift uncertainties Figure 4 shows the 
increase in error bars when marginalizing over both IA and photo-z 
uncertainties. We model the latter assuming that a Gaussian photo-z 
probability distribution function is obtained from each galaxy. Un¬ 
certainties in the overall redshift distributions of each tomographic 
bin are then parameterized as variations in the mean (“bias”) and 
variations in the width of said Gaussian “cr”. We further assume 
Gaussian priors on the bias and cr parameters, i.e. we consider an 
optimistic photo-z scenario with Abias = 0.002, Acr = 0.003, and a 
pessimistic scenario with Abias = 0.005, Act = 0.006. 

The red/dashed contours in Fig. 4 show results when marginal¬ 
izing over photo-z uncertainties only. We see a substantial degrada¬ 
tion of constraints on fl ra , <x 8 , and to a lesser extent on Wo- Ac¬ 
counting for IA and photo-z jointly also substantially degrades w a 


© 0000 RAS, MNRAS 000. 000-000 





































The impact of intrinsic alignment on current and future cosmic shear surveys 9 


and n s . In the presence of IA the two levels of photo-z errors we 
consider hardly change the contour size. When looking at the nui¬ 
sance parameter panel of Fig. 4 we see that although the photo-z 
parameters (bias ph , and cr ph ) are obviously better constrained in the 
optimistic photo-z prior scenario, the additional information gets 
largely absorbed in the luminosity function parameters. 

We emphasize that the Gaussian photo-z model is optimistic 
in the sense that it neglects catastrophic outliers, which have been 
shown to severely degrade dark energy constraints (Bernstein & 
Huterer 2010). We postpone a more quantitive analysis of the de¬ 
generacies of IA, photo-z, and cosmological parameters to future 
work. 

Varying the IA model We now extend our analysis to all IA sce¬ 
narios described in Sect. 2.2. Specifically, we create IA contami¬ 
nated data vectors for the NLA Cosmic emulator, the Freeze-in, the 
linear, and the tidal alignment model, and then analyze these data 
vectors assuming the NLA Halofit model in the marginalization. 
This analysis mimics the realistic scenario where the analyst has a 
broad idea, but not an exact model of the IA contamination in the 
data. 

Figure 5 shows the 68% likelihood contours for the aforemen¬ 
tioned four IA scenarios. The NLA Cosmic Emulator (black/solid) 
and the tidal alignment (green/long-dashed) scenarios are well cap¬ 
tured by the NLA Halofit model. The linear power spectrum sce¬ 
nario (red/dashed) shows the largest offset from the underlying 
fiducial model, the freeze-in scenario’s offset (blue/dot-dashed) is 
smaller but noticeable. Even for the worst case, i.e. the linear IA 
contamination mitigated using a NLA model, the fiducial cosmo¬ 
logical model is within the 68% likelihood contours. As our like¬ 
lihood analysis assumes / max = 5000, where differences between 
linear and non-linear power spectra are significant, we conclude 
that when IA is fully described by some tidal alignment scenario, 
an exact IA model is not vital for removing biases in cosmological 
parameters. We will consider different surveys and the impact of 
the galaxy luminosity function in the following subsection. 

5.3 Comparison of LSST, DES, Euclid, WFIRST 

In the previous sections we examined the impact of IA contami¬ 
nation on LSST WL constraints; we now extend our analysis to 
WFIRST, Euclid, and DES. For the comparison of these surveys it 
is important to note that the significance of IA biases, and hence 
the importance of IA as a contaminant, is affected by survey depth 
and survey width: The amplitude of the IA contamination relative 
to the cosmic shear signal is highest at low redshift and for more lu¬ 
minous source galaxies, hence the amplitude of the parameter bias 
decreases with survey depth. The significance of the contamination 
is given by the ratio of the amplitude of the parameter bias and the 
survey’s constraining power on these parameters, which increases 
with survey depth and width. 

Figure 6 shows the marginalized one-dimensional 68% er¬ 
ror bars on cosmological parameters for the deep surveys LSST 
(m r = 27) and WFIRST (m r = 28). Figure 7 shows the same 
analysis for the shallower surveys Euclid (m r = 24.5) and DES 
(m, = 24). Each plot comprises results of 22 different likelihood 
analyses; the presentation of the results for different surveys follow 
the same structure: 

We first show the statistical error bars (black) obtained from 
using an uncontaminated DM data vector as an input and switching 
off the CosmoLike IA module in the simulated likelihood analysis. 
In blue we show a similar analysis but reducing the number density 


of source galaxies by 20% (denoted as “BL”), which is a simplified 
implementation of identifying the severely IA contaminated (red) 
fraction of galaxies and excluding these galaxies from the analysis 
(c.f. Fig. 2 for the expected red fraction of different surveys). 

We note that the reduction in source density only causes a 
minor loss of cosmological constraining power (also see Fig. 8). 
While this result might be counter-intuitive given the apparent im¬ 
portance of shape noise in many forecasts, it can be explained by 
the use of non-Gaussain covariance matrices. For Gaussian covari¬ 
ances only the cosmic variance term is independent of rc ga i, whereas 
a larger n ga i substantially reduces the noise and mixed term, which 
again substantially increases the information on small physical 
scales. Non-Gaussian covariances are dominated by terms that stem 
from higher-order moments of the density field and especially the 
Halo Sample Variance term, all of which are unaffected by n gal . This 
result emphasizes the prospect of selecting well-understood and 
characterized galaxy subsamples for the analysis, thereby avoiding 
potentially severe systematics while minimally degrading cosmo¬ 
logical constraints. 

After establishing the statistical error budget as a baseline we 
consider five IA scenarios (see Sect. 2.2 for details) to contaminate 
the data vector. In order of their appearance in Figs. 6 and 7 the 
considered models are: 

• Nonlinear linear alignment using Halofit P s (orange) 

• Nonlinear linear alignment using Coyote Universe P s (red) 

• Linear alignment using linear power spectrum (transfer func¬ 
tion from Eisenstein & Hu 1998) (green) 

• Freeze-in model using the geometric mean of linear power 
spectrum and Halofit Ps (blue) 

• Full tidal alignment model using the Halofit P s (dark blue) 

For each model we run four different analyses, two for the 
GAMA LF and two for the DEEP2 LF. Of these two the first anal¬ 
ysis always shows the impact of the IA contaminated data vector in 
the absence of any mitigation and the second analysis shows con¬ 
straints when using the NLA Halofit model in the CosmoLike IA 
mitigation module. The marginalized likelihood is obtained by in¬ 
tegrating over a ten-dimensional IA nuisance parameter space as 
described in Eq. (30). 

In addition to Figs. 6 and 7, we evaluate the magnitude of 
the bias in cosmological parameters from all likelihood analyses 
using a similar metric as in Eifler et al. (2014a). The bias is de¬ 
fined as the difference between the best-fit value of the parameter 
emerging from the likelihood analysis and the fiducial value used 
to generate the spectra. We then compare this bias to the size of 
the seven-dimensional cosmological likelihood volume at different 
confidence levels and asses whether our IA mitigation strategy is 
successful. 

If we considered only one cosmological parameter the relevant 
quantity that characterizes bias with respect to likelihood volume 
would be (p bestfit _ p fld ) 2 /cr 2 . For our seven parameter case, this 
concept generalizes to 

At 2 = (Pc d - p' A ' beStflt )' C-‘ (p“ - pIA.bestflt)^ (32) 

where the covariance matrix is determined via 

<£ = WTT Z (< p ‘> ~ p *) (< p ‘> - p 0 ’ (33) 

k=0 

with (pj.) indicating the mean of the z'-th cosmological parameter 
(i,j e [1,7]), and k e [1,1V] being the index running over all steps 
in the MCMC chain. Assuming a^ 2 distribution with seven degrees 
of freedom, we find the critical tsf 1 values that correspond to 68%, 
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Figure 5. Top: WL constraints (68% confidence region) from LSST when assuming different models in the data vector, but assuming the fiducial IA NLA 
Halofit model in the marginalization. The black/solid lines show results from a data vector contaminated by the NLA Coyote Universe model, which is very 
close to the fiducial NLA Halofit scenario. Red/dashed corresponds to the linear power spectrum model, blue/dot-dashed to the freeze-in model, green/long- 
dashed to the tidal alignment model. All models assume the GAMA luminosity function. Bottom: The marginalized one-dimensional posterior probabilities 
for the ten nuisance parameters describing the IA uncertainty. 


95%, and 99% confidence regions are 8.14, 14.07, and 18.48, re¬ 
spectively. We summarize corresponding values in Table 5 and find 
a qualitative similar behavior as in the Figs. 6 and 7. 

It is important to note that Ay 2 and the metric employed in 
Figs. 6 and 7 quantify the success of IA mitigation differently. The 
latter accounts for the non-Gaussianity of the posterior distribu¬ 
tion in parameter space, but projects the posterior probability onto 
a specific basis (i.e. the cosmological parameters). It is fair to as¬ 
sume that if all projections are free of residual biases, the overall 
bias for the seven-dimensional cosmological parameter space is un¬ 
der control. The A^ 2 values in Table 5 give an overall estimate of 
the residual bias in combination with the multi-dimensional error 


in parameter space, however under the explicit approximation that 
the posterior probability is distributed as a multivariate Gaussian. 

IA impact For a given survey, the non-linear IA models (NLA 
Halofit, NLA Coyote, TA) in combination with the DEEP2 lumi¬ 
nosity function give rise to the strongest IA contamination, and 
hence cause the largest biases in cosmological parameters when 
no mitigation is applied. 

As a consequence of small statistical uncertainties and rela¬ 
tively shallow observations, Euclid is most severely affected by IA. 
Corresponding A^ 2 values in Table 5 range from 111 for the NLA 
Halofit to 91 for the TA scenario. While LSST has similar con¬ 
straining power as Euclid, the relative contribution of the IA con- 
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Figure 6. One-dimensional 68% error bars on cosmological parameters for the deep surveys LSST and WFIRST. The abbreviations that reference the different 
models read “DM” :=Dark matter, “BL” := DM with Blue Galaxies only, “H” := NLA Halofit, “C”:=NLA Coyote Universe, “L”:= LA linear power spectrum, 
“F” : _ f reeze _j n “T”;= tidal alignment. For each model we run four different analyses, considering impact (”i”) of untreated IA contamination and marginalized 
(“m”) results using the NLA Halofit model to account for the IA uncertainty. For all IA scenarios we consider both GAMA and DEEP2 luminosity functions, 
denoted as “G” and “D”, respectively. 


tamination to the shear signal is reduced due to the deeper redshift 
distribution, and the corresponding parameter biases are smaller. 
For LSST the Ax 2 values for these IA scenarios range from 71 to 
57. 

DES has substantially less constraining power than either Eu¬ 
clid or LSST, but it is also significantly shallower than LSST, and 
even slightly shallower than Euclid; as a consequence, despite the 
larger error bars, parameter biases due to unmitigated IA are still 
considerable and Af 1 values for the non-linear IA models range 
from 58 to 54. The HLS component of the WFIRST survey with an 


area of 2200 deg 2 has higher constraining power than DES, but less 
than Euclid and LSST; due to the depth of the survey, A^ 2 values 
for these three IA scenarios range from 37 to 16, which still consti¬ 
tutes a significant bias that must be addressed. 

Compared to the non-linear IA models, the IA contamination 
from the linear alignment and freeze-in model is much smaller (at 
large I), and the parameter biases are less severe. For WFIRST these 
biases are below the 1-cr threshold before any mitigation is applied, 
i.e. A^ 2 = 4.1 and A^ 2 = 2.3 for freeze-in and linear alignment, 
respectively. 
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Figure 7. Same as Fig. 6, but for the shallow surveys DES and Euclid. 


IA mitigation The CosmoLike IA mitigation module (NLA 
Halofit model with 10 nuisance parameters, c.f. Sect. 5.1) success¬ 
fully removes IA biases for the NLA Halofit, NLA Coyote, and the 
TA scenarios. While this is expected for the NLA Halofit scenario 
as input (recall that we use this exact scenario in our mitigation 
module), it is encouraging to see that the technique is sufficiently 
robust to mitigate the other two scenarios with different prescrip¬ 
tions for the non-linear regime as well. In particular, biases caused 
by the full tidal alignment model, which as Blazek, Vlah & Seljak 
(2015) have shown is in excellent agreement with observations, are 
fully removed by our mitigation scheme. The corresponding A y 1 
values for this model assuming the DEEP2 LF and Euclid survey 
drop from 91 to 1.1 after mitigation, indicating that after IA mitiga¬ 


tion the underlying cosmology can be recovered without significant 
biases. 

For the linear alignment and freeze-in models mitigation is not 
as straightforward: These models predict angular IA power spectra 
which differ substantially in terms of shape in / and redshift scaling 
compared to the NLA Halofit template assumed in our mitigation 
strategy. As a consequence the success of the IA mitigation de¬ 
pends on how well the nuisance parameters can compensate these 
differences. In particular, if differences between these IA models 
are more degenerate with cosmological parameters than with nui¬ 
sance parameters the IA mitigation will not be successful. 

For the freeze-in model IA mitigation using the NLA Halofit 
works sufficiently well for the LSST, WFIRST and DES survey pa¬ 
rameters. For Euclid however, the residual bias may be severe de- 
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pending on the underlying luminosity function: we find A^ 2 = 19 
and A^ 2 = 7.4 for the DEEP2 LF and GAMA LF parameters re¬ 
spectively. This indicates that the simple mitigation template pro¬ 
posed in the analysis may not be sufficient for Euclid if the IA con¬ 
tamination were similar to the freeze-in model. This is clearly a 
consequence of the fact that we use use the NLA Halofit model in 
our IA mitigation strategy independent of which scenario was as¬ 
sumed in the data vector. Even though the freeze-in model is based 
on a combination of linear and nonlinear power spectrum, its an¬ 
gular shape and redshift dependence are sufficiently different from 
the NLA Halofit that ten (amplitude related) nuisance parameters 
cannot compensate these differences. 

This problem becomes more prominent for linear alignment 
model as the underlying IA scenario. For Euclid, the residual bias 
after IA mitigation is significant for the GAMA LF model, and even 
more severe for the DEEP2 LF (corresponding values are = 11 
and A^ 2 = 29, respectively). Even for LSST, the DEEP2 LF case 
after IA mitigation (A^ 2 = 11) now exceeds the 1-cr threshold 
(A^ 2 = 8.14). The difference between linear and NLA IA contami¬ 
nation, and thus the failure of our mitigation strategy, is aggravated 
by the optimistic choice of scale ( max = 5000. For both, Euclid and 
LSST, we run additional chains with a reduced / max = 2000; the 
corresponding tsf 2 values drop to 13 (from 29) and to 6 (from 11) 
for Euclid and LSST, respectively. This allows us to claim moderate 
success for LSST, Euclid however is still severely affected. 

For WFIRST, the proposed IA mitigation is successful; there 
was no large bias to begin with and mitigation shows a slight im¬ 
provement, primarily due to a slightly increased error bar after 
marginalizing over the nuisance parameters. A similar statement 
may hold for DES, however the bias in case of the DEEP2 LF ap¬ 
proaches the 1-cr threshold. We discuss the implications of these 
findings on IA mitigation strategies for future surveys in Sect. 6. 

In addition to the bias removal we consider the information 
loss as a major criteria for a successful mitigation scheme. We 
quantify this information loss “IL” as the ratio of geometric mean 
computed from the 1-D cosmological parameter error bars of the 
considered scenario and the geometric mean statistical error. 


IL = 


n, 


n,of 


1/7 


(34) 


Table 6 shows the corresponding increase in error bars for five dif¬ 
ferent scenarios, i.e. when removing 20% of the galaxies, the NLA 
Coyote scenario, and the full tidal alignment scenario. The loss in 
constraining power when removing one fifth of the galaxies is min¬ 
imal across all surveys; it ranges from 4-10%. These numbers are 
higher when mitigating IA with 10 nuisance parameters and pri¬ 
ors as described in Tab. 4; the increase in error bars ranges from 
approximately 30% to 55% depending on survey and luminosity 
function. 

For the interpretation of these numbers it is important to note 
that this analysis employed broader priors range than current obser¬ 
vational limits for the NLA amplitude (e.g. Joachimi et al. 2011) 
and LF parameters, in order to enable the comparison of differ¬ 
ent IA models, such as the mitigation on linear alignments with 
a halofit NLA template, which requires normalization parameters 
outside the observed prior range for the NLA model. Hence the 
losses in constraining power reported in Table 6 are conservative 
estimates, the use of observed priors for the NLA amplitude and 
LF measurements from the data set will improve the performance 
of this mitigation strategy considerably. 


Table 5. The Ay 2 distance (see Eq. 32) between best fit and fiducial param¬ 
eter point. Cases where IA mitigation fails are indicated as red. 


Scenario 

DES 

Euclid 

WFIRST 

LSST 

DM 


0.39 

0.15 

0.35 

0.094 

BL 


0.34 

0.21 

0.38 

0.16 

HGi 


22 

65 

6.6 

61 

HGm 


0.56 

0.36 

0.61 

0.42 

HDi 


55 

111 

18 

73 

HDm 


0.59 

0.27 

0.75 

0.29 

CGi 


22 

63 

5.9 

55 

CGm 


0.82 

0.44 

0.79 

0.36 

CDi 


58 

99 

16 

71 

CDm 


0.54 

0.42 

0.73 

0.42 

LGi 


3.4 

15 

0.94 

5.8 

LGm 


2.1 

11 

0.55 

2.5 

LDi 


11 

33 

2.3 

15 

LDi (/max = 

2000) 

- 

39 

- 

19 

LDm 


5.9 

29 

1.4 

11 

LDm (/max 

= 2000) 

- 

13 

- 

6 

FGi 


5.8 

32 

1.4 

14 

FGm 


1.8 

7.4 

0.82 

1.8 

FDi 


19 

68 

4.1 

33 

FDm 


3.3 

19 

0.91 

3.8 

TGi 


30 

87 

13 

64 

TGm 


1.5 

0.8 

0.55 

0.69 

TDi 


54 

91 

37 

57 

TDm 


1 

1.1 

0.99 

0.77 


Table 6. The loss in constraining power with respect to the statistical errors 
for five scenarios (see text for details). 


Scenario 

DES 

Euclid 

WFIRST 

LSST 

BL 

1.05 

1.07 

1.04 

1.1 

CGm 

1.42 

1.44 

1.4 

1.39 

CDm 

1.35 

1.48 

1.49 

1.52 

TGm 

1.42 

1.53 

1.37 

1.49 

TDm 

1.33 

1.55 

1.39 

1.51 


5.4 Alternative IA mitigation strategies 

In this section we investigate alternative IA mitigation strate¬ 
gies for Euclid, since we have established that Euclid is most 
severely impacted by IA. Figure 8 shows the Euclid statistical er¬ 
rors ( black/solid) in comparison to the bias that occurs when an¬ 
alyzing the fiducial Euclid IA data vector without any mitigation 
{red/dashed). 

Under the assumptions that (1) only red galaxies are intrinsi¬ 
cally aligned and that (2) we can precisely distinguish between red 
and blue galaxies, it is an obvious mitigation strategy to remove the 
red fraction. The blue/dot-dashed contours correspond to removing 
20% of the Euclid galaxies, which results in a number density of 
16/arcmin 2 . The constraining power is hardly diminished, which as 
mentioned before can be explained by the fact that the covariance is 
dominated by non-Gaussian terms rather than the shape noise term. 
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Figure 8. Impact of IA (NLA Halofit, GAMA LF) for the Euclid mission ( red/dashed ) compared to the Euclid statistical errors ( black/solid ). The blue/dot- 
dashed contours correspond to removing 20% of the Euclid galaxies, which mimics an IA mitigation through conservatively removing the contaminated, red 
galaxy fraction. The green/long-dashed contours show the impact of IA (fiducial model,/ re( i = 20%), but with IA being completely removed for z < 1, which 
mimics an analysis that includes an overlapping low-z spectroscopic survey. 


Obviously the assumptions (1) and (2) are optimistic. Tidal 
torquing mechanisms cause higher-order IA effects in blue galax¬ 
ies, and on sufficiently large scales, a contribution linear in the tidal 
field may dominate for these galaxies as well. Given the current 
observational upper limits on IA of blue galaxies, the effect is non- 
negligible and can potentially cause severe biases in cosmological 
parameters. Targeted spectroscopic information is needed to im¬ 
prove our understanding of IA of these galaxies. 

The green/long-dashed contours correspond to an analysis 
where we assume an IA contaminated (red galaxies only) Euclid 
data vector (NLA Halofit, GAMA LF), but the IA signal is set to 
zero for z < I. We find that the severe bias seen in red/dashed con¬ 
tours has almost completely vanished, indicating that at high red- 
shifts the shear signal completely dominates the IA contamination. 

This result motivates a low-z spectroscopic survey (or multi¬ 
band photometric survey as proposed in Dore et al. (2014)) that 
overlaps with the corresponding imaging survey to determine the 
red galaxy IA scale dependence and amplitude over the relevant 
redshift range. If a negligible level of IA from blue galaxies (con¬ 
sistent with current non-detections) can be confirmed with tar¬ 
geted spectroscopic observations, a wide, shallow (z < 1) spectro¬ 


scopic survey may be sufficient for the calibration of IA mitigation 
schemes. 


6 DISCUSSION AND CONCLUSIONS 

In this paper we present forecasts on the impact of intrinsic align¬ 
ments on future surveys and explore the performance of an IA mit¬ 
igation strategy that uses the NLA model as a template, employing 
10-12 nuisance parameters to account for luminosity and redshift 
scaling of the IA amplitude and the color and redshift distribution 
of source galaxies. 

Compared to previous studies of the impact of IA on cosmic 
shear (e.g., Kitching, Taylor & Heavens 2008; Joachimi & Bridle 
2010; Kirk, Bridle & Schneider 2010; Kirk et al. 2012, 2015), this 
analysis uses more detailed and realistic models for the IA contam¬ 
ination, the expected statistical noise, and forecasting methods: 

• In this study, modeling of the IA amplitude accounts for 
survey specific red fraction and luminosity distribution of source 
galaxies, while previous studies (e.g., Kirk et al. 2015) used the 
amplitude normalization of the bright SuperCOSMOS galaxy sam¬ 
ple (Hirata & Seljak 2004; Bridle & King 2007) to represent the 
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alignment of all source galaxies (but see Joachimi et al. 2015, for 
a recent forecast including luminosity dependent normalization). 
We note that accounting for the luminosity dependence of the IA 
amplitude reduces the amplitude of parameter biases considerably. 

• Previous forecasts relied on Gaussian covariances while we 
include higher-order moments of the density field in our error bars. 
In particular the inclusion of Halo Sample Variance causes a signif¬ 
icant correlation of power spectrum measurement at large l, which 
increases the statistical uncertainty and hence decreases the signif¬ 
icance of biases. 

• Previous analyses used the Fisher forecasting formalism, 
while this study uses simulated likelihood analyses using MCMC. 
Fisher analyses inherently assume that the surface of the posterior 
probability is described by a Gaussian, and the systematic effects 
can be linearized in the nuisance parameters to calculate biases. 
While this is an accurate description at the best-fit parameter point, 
it becomes much less accurate for the outer regions of the parameter 
space. While it is not clear whether Fisher analyses over- or under¬ 
estimate the bias (it can go either way), it is clear that an MCMC 
is more accurate and may show very different results, especially in 
high-dimensional, degenerate parameter spaces. 

• We consider a large range of IA scenarios to encompass IA 
modeling uncertainties. We also consider the recently developed 
full tidal alignment model Blazek, Vlah & Seljak (2015), which has 
been found to be in excellent agreement with IA measurements. 

Based on these improved forecasting capabilities we find the 
unmitigated bias from IA on cosmological parameters to be less 
severe than most previous studies. The significance of these biases 
is highest for shallow surveys (as the relative IA contamination is 
highest at low redshift), and increases with survey area as the sur¬ 
vey’s constraining power increases. 

As described in Sect. 4.1, we average the observed redshift 
and luminosity scaling on the IA amplitude from Joachimi et al. 
(2011) over the source galaxy LF. Due to the current lack of repre¬ 
sentative luminosity function measurements to the depth required 
for deep future surveys, we extrapolate two moderate-redshift mea¬ 
surements of the LF and its evolution. The difference between these 
two models is significant, and modeling of the galaxy LF should be 
considered a key uncertainty in future IA forecasts. 

We propose an IA mitigation strategy that is based on the NLA 
Halofit template for the IA power spectrum shape and uses 10 nui¬ 
sance parameters to model the amplitude as a function of redshift 
and source galaxy distribution. We find this mitigation strategy to 
be sufficient to remove the parameter bias for the non-linear IA 
models considered in this study (NLA Halofit, NLA Coyote, full 
tidal alignment), and for all surveys (DES, Euclid, LSST, WFIRST) 
considered. This indicates that for the currently planned stage IV 
surveys IA mitigation does not require an exact model for IA in the 
non-linear regime; an approximate model will be sufficient. 

The proposed mitigation strategy breaks down if the differ¬ 
ence between the IA scenario contaminating the data vector and 
the NLA Halofit template assumed in the mitigation becomes too 
large, as demonstrated by the example of mitigating a linear align¬ 
ment contamination with a NLA Halofit template. However, recent 
comparisons of observational data and IA models (e.g. Joachimi 
et al. 2011; Blazek, Vlah & Seljak 2015) strongly prefer NLA and 
full tidal alignment models over the linear alignment model, hence 
this specific failure is likely less of a concern for practical purposes. 
However, in order to test the fixed shape assumption, one can in¬ 
clude a simple parameterization of the IA power spectrum shape, 


e.g. Pse(k, z) = P^(k, z)P\lf\k, z) with /3 as an additional nuisance 
parameter, in future implementations of this mitigation scheme. 

As implemented here, this mitigation strategy reduces a sur¬ 
vey’s constraining power on individual cosmology parameters by 
several tens of percents (c.f. Table 6). We emphasize that these 
numbers are overly pessimistic due to the priors on IA and LF pa¬ 
rameters employed in this analysis (see Table 4). These are sev¬ 
eral times broader than current observational constraints, which 
was required due to the broad range of IA power spectrum models 
presented here. It it unclear how well these parameterization will 
extend to faint luminosities and high redshifts, hence measuring 
the IA parameters as demonstrated in Mandelbaum et al. (2011); 
Joachimi et al. (2011), and the LF of the source galaxy sample is 
the most promising way forward to further enhance IA mitigation 
schemes. The wide prior ranges on IA and LF parameters are also 
the reason why our joint analysis of IA and Gaussian photo-z uncer¬ 
tainties is dominated by the former. We defer an in-depth analysis 
of the degeneracies between IA self-calibration with stringent pri¬ 
ors and more realistic photometric redshift uncertainties to future 
work. 

In general, the success of a mitigation scheme that marginal¬ 
izes over the parameters of a specific model for a systematic ef¬ 
fect will depend on whether the true signal from this effect is more 
degenerate with the chosen nuisance parameters (success) or the 
cosmological parameters (failure). If multiple systematics mod¬ 
els are viable, one might consider a hybrid model, e.g. allowing 
the analysis code to use a weighted combination of two models, 
A x Model] + B x ModeL, and marginalizing over A,B in addi¬ 
tion to the nuisance parameters within Model] and ModeL. In the 
context of IA mitigation, this will be particularly relevant for more 
detailed modeling and mitigation of blue galaxy alignments. 

If the upper limits on blue IA continue to decrease with future 
observations, and if red and blue galaxies can be identified with suf¬ 
ficient accuracy, our analysis suggests removing red galaxies from 
the cosmic shear analysis. One interesting results of this paper is 
that a removal of 20% of the galaxies only cause a loss in infor¬ 
mation of 4-10% depending on the survey considered. Choosing a 
well-understood galaxy (sub)sample for future cosmic shear anal¬ 
yses, instead of maximizing the galaxy number density at all costs, 
may be a favorable strategy. 

Any IA mitigation technique will be further enhanced by ex¬ 
tending the cosmic shear analysis to a joint analysis of cosmic shear 
and other observables required to further self-calibrate the IA am¬ 
plitude from the same data set, i.e. low-redshift measurements of 
the (photometric) galaxy position-shape correlation, and cluster¬ 
ing correlation function measurements for the same galaxy sample 
(to self-calibrate galaxy bias parameters required for modeling the 
galaxy position-shape correlation). However, such a joint analysis 
including non-Gaussian (cross-) covariances for all observables is 
beyond the scope of this paper. 
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